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ABSTRACT 


The muon g-2 experiment at Brookhaven National Laboratory mea¬ 
sures the anomalous magnetic moment of the muon, a /t , very precisely. 
This measurement tests the Standard Model. The analysis for the data 
collected in 2000 (a /i + run) is completed and the accuracy on is 
improved to 0.7 ppm, including statistical and systematic errors. The 
data analysis was performed blindly between the precession frequency 
and the field analysis groups in order to prevent a possible bias in the 
result. The result is a A( (exp) = 11 659 204(7)(5) x 10 -10 (0.7 ppm). 
This paper features a detailed description of one of the four analyses 
and an update of the theory. 



1 Introduction 


The gyromagnetic ratio (g-factor) of a particle is defined as the ratio of its mag¬ 
netic moment to its spin angular momentum. For a point-like spin-1/2 particle, 
the g value is predicted by the Dirac equation to be equal to 2. On the other 
hand, from experiments on the hyperfine structure of hydrogen in the late 1940’s, 
it was found that the g value of an electron was not exactly 2 (« 2.002). In 
the first order this deviation is clue to the creation and absorption of a virtual 
photon by the particle. This process is described by quantum electrodynamics 
(QED) and the electron g-2 experiments were precise tests of QED. The Stan¬ 
dard Model (SM) calculations for the anomalous magnetic moment of the muon 
includes also the small contributions from hadronic and weak interactions, in ad¬ 
dition to electromagnetic (QED) interactions. The QED contribution to the SM 
is the largest, however it is the most precisely known. On the other hand, the 
largest contribution to the uncertainty comes from the hadronic interactions. 

2 Experimental Method and Apparatus 

When a positive muon decays into a positron and two neutrinos, parity is max¬ 
imally violated so that the average positron momentum is along the muon spin 
direction (in the muon rest frame). Since observing the decay energy in the labora¬ 
tory frame is equivalent to observing the decay angle in the center-of-mass system, 
the muon spin direction can be measured by counting the decay positrons above a 
certain energy threshold. As a result, the positron counting rate is modulated with 
the muon spin precession. This is the main concept of the muon g-2 experiments. 1 

In our experiment, polarized muons are injected into a superconducting storage 
ring. The ring provides a very uniform static dipole magnetic field, B. If the 
g-factor of a muon were exactly equal to two, the spins and the momenta of the 
muons would stay parallel during the storage. This would mean the cyclotron 
frequency, u c , of the muon is equal to the spin precession frequency, ca a . However, 
due to the anomaly on the magnetic moment, the spin precesses faster than the 
momentum, uj s = uy(l + ayy). The muon spin precession angular frequency 
relative to its momentum in the presence of a vertical focusing electric field E and 
magnetic field is given by: 
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where a p is the muon anomalous magnetic moment and 7 is the relativistic Lorentz 
factor. 


The influence of the electric held on uj a is eliminated by using muons with the 
“magic” momentum 7 M = 29.3, for which P = 3.09 GeV/c. The measurement of 
the angular precession frequency u: a and static held B measured in terms of the 
proton NMR frequency, uj p , determines a p as follows; 
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where A = = 3.183 345 39(10 ). 2 


The counting rate N(t) of decay positrons with energies above an energy threshold 
E is modulated by the spin precession of the muon, ideally leading to 

N{t) = N 0 (E) e~ t/T » (1 + ALcos [u a t + </>]) (3) 

where N 0 is the normalization, r /t is the time-dilated muon lifetime, A is the 
energy-dependent asymmetry and (j) is the phase. The angular precession fre¬ 
quency u a is determined from a ht to the experimental data. Figure 1 shows the 
g-2 precession frequency data collected in 2000. The error bars are in blue. Figure 
2 shows the comparison of the energy spectrum at the peak and the trough of the 
g-2 cycles. For a certain energy threshold, the difference between two energy spec¬ 
tra is reflected in the asymmetry parameter we measure. The energy dependence 
of the asymmetry can be easily seen. To minimize the statistical error on the pre¬ 
cession frequency, the energy threshold is chosen between 1.8-2.0 GeV. For this 
threshold, the number of events times the square of the asymmetry is maximum. 
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Fig. 1 : g-2 time spectrum. The pe¬ 
riod of the g-2 is 4.36 (is and the dilated 
lifetime is 64.4 /is. 


Fig. 2 : A comparison of the posit¬ 
ron energy spectrum obtained at the 
peak and at the trough of g-2. 


The BNL g-2 experiment uses a polarized muon beam produced by the Alter¬ 
nating Gradient Synchrotron (AGS), which is the brightest proton source in its 
energy class in the world. Protons with an energy of 24 GeV hit a nickel target; 
pions from this reaction are directed into a pion decay channel where they decay 
into muons. Muons are naturally polarized when a small forward momentum bite 
is taken. 

Our storage ring is a single superconducting C-shape magnet, supplying a 1.45 
T uniform field. 3 Beam muons are injected into the storage ring through a su¬ 
perconducting inflector magnet, which locally cancels the storage ring field and 
delivers muons approximately parallel to the central orbit. 4 

The BNL g-2 experiment started first data taking in 1997 using pion injection 
into the storage ring. 5 This was similar to what was done in the last CERN ex¬ 
periment. 1 Since the pions fill most of the phase space, there are always a few 
daughter muons from pion decay which would be captured into stable equilibrium 
orbits in the storage ring. The disadvantage of pion injection, was the high back¬ 
ground level due to the pions themselves. From 1998 on, our experiment used 
muon injection. In order to put the muons onto their equilibrium orbit, a fast 








magnetic kicker was used. The background level was reduced dramatically and 
the number of events was substantially increased with muon injection compared 
to pion injection. The kicker is located at 90° with respect to the inflector and 
consists of three sets of parallel plates each 1.7 nr long, carrying 5200 A for a very 
short time (basewidth 400 ns). 6 

Electrostatic quadruples 7 are used to vertically confine the muons. Four quadrupole 
assemblies are located in the ring, each pulsed with ±24 kV and covering 43% of 
the azimuth in the storage ring. 

Inside the ring, the decay positrons are observed for approximately 600 /xs by 24 
electromagnetic shower calorimeters 8 consisting of scintillating fibers embedded in 
lead. All positrons above a certain threshold are digitized individually with 400 
MHz waveform digitizers and the digitized waveforms are stored for the analysis. 

3 uj a Analysis 

The g-2 frequency u a is determined by fitting the time spectrum of positrons after 
data selection. There were four independent analyses of the precession frequency 
from the 2000 data. However, only one of them 9 is going to be described here. 

When statistics are high, the influence of small effects becomes observable in 
the data. Deviations in the y 2 and the parameter stability when experimental 
conditions are varied are the signature of the size of these effects. The 2000 data 
amount to four times the number of events recorded in 1999. Positron pilcup, 
coherent betatron oscillations (CBO) and muon losses were already observed in 
1999 and accounted for in the fitting function. 10 The mismatch between the 
inflector and storage ring acceptances is one of the main sources of CBO. CBO is 
also caused by a non-ideal kick to muons when they first enter the ring. Therefore, 
they are not captured in their equilibrium orbit, but rather oscillate about it. As 
a result of these oscillations, the beam gets closer and further away from the 
detectors resulting in a modulation of the counting rate. CBO can be easily seen 
in the residuals constructed from the difference between the data and the ideal 
5-parameter fit (Eq.3 and Fig. 3). The Fourier transform of the residuals is shown 
in Figure 4. 
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Fig. 3 : Time spectrum of the residuals. Fig. 4 : Fourier spectrum of the residuals. 


Modulations due to CBO are parameterized as: 

Fcboif) = 1 + A cbo e ~ t,T ^° cos [2ir f cbo t + 0 c6o ] (4) 

where A cbo is the amplitude of the modulation, r cbo is the coherence time of the 
damping, and (f) cbo is the phase. The frequency f cbo is fixed to the frequency 
determined by the Fourier analysis (465.7±0.1 kHz) of the residuals. With the 
larger data set from 2000 , additionally the g -2 asymmetry (H) and the phase (0) 
are also seen to be modulated with the CBO. Therefore, the fitting function (Eq. 
3) was modified as follows : 

N(t) = N 0 F cbo {t) e~ t/r » (1 + A(t) cos [u a t + 0(f)]) (5) 

where 

A(t) = A ( 1 + A A e ~ t/Tcb ° cos [iv cbo t + cp A ] ) ( 6 ) 

and 

0(f) = 0 + e~ t/Tcbo cos [ujchot + 0<^>] • (7) 

The amplitudes A cbo , A a and A $ are small (« 1%, ~ 0.1% and ~ 1 mrad, respec¬ 
tively) and are consistent with Monte Carlo simulation. 

To avoid beam resonances, the weak-focusing index n was set at 0.137, half way 
between two neighboring resonances (Figure 5). Running at a higher n value 


































was not possible because the high voltage on the quads was limited to avoid 
breakdown. Running at a lower n value leads to storing less beam. Running at 
n = 0.137 made the CBO frequency close to two times the g-2 frequency, which 
was a difficulty for the analysis, especially for the fit. 


Pilcup, the overlap of positron signals 
in time, leads to a misidentihcation of 
positrons, their arrival time and their en¬ 
ergy. Pilcup events have half the lifetime 
of the muons. The size of the effect is 
larger at early times compared to late 
times. The way the data was stored per¬ 
mitted us to resolve pilcup using the data 
itself, ft was preferred to subtract the 
pilcup events from the data because in¬ 
cluding the pileup related terms into the 

Fig. 5 : Beam tune and n value. The 
lines represent the resonances. 

fitting function caused cross-talk between the fit parameters and an increase by a 
factor of two in the uncertainty of u a . Pileup pulses were reconstructed artificially 
from the data itself in the following way. Every positron pulse above a relatively 
high threshold is digitized at 400 MHz for 80 ns, which is called a “WFD island”. 
Positron pulses are fitted with a pulse-finding algorithm to determine energy and 
time information. This pulse-finding algorithm can resolve only events that are 
more than 3 ns apart. After the main triggering pulse, there may be a second 
pulse on the WFD island, which carries the necessary energy and time informa¬ 
tion for artificial pileup construction. When there is a main triggering pulse, we 
looked for a secondary pulse on the same WFD island within a time window off¬ 
set from the trigger. If there is a pulse, the main and the secondary pulses are 
added properly and the time is assigned from the energy-weighted time of the pair 
(Figure 6). These constructed pileup events are then subtracted from the data 
to obtain a pileup-free time spectrum. Figure 7 shows the fit to constructed pileup. 











Pileup Candidates Resolution Time = 2 , 



Constructed Doubles 


E=E1+E2 



<t>=<tl,t2-Offset Time> 


D 


Fig. 6 : Schematic view of the pileup 
construction. 



Fig. 7 : Fit to constructed pileup. 


One of the tests of checking the pileup subtraction quality is to look at the average 
energy versus time. The average energy of the detected positrons, for a given g-2 
period, with and without the pileup subtraction, is shown in Figure 8. The effect 
of pileup can be seen clearly in the upper curve. After the pileup subtraction the 
average energy versus time is flat (lower curve). 

The next step in the analysis was to include muon losses. The muon losses were 
determined two different ways. The pileup subtracted data were fitted at later 
times (~ 300/is) to an 8-parameter functional form (Ideal+CBO) and extrapo¬ 
lated to earlier times. Only the CBO related parameters were determined at 50 /is 
since the lifetime of the CBO is ~ 110/rs. The ratio of the data to the extrapolated 
fitting function shows the effect of the lost muons. The second method is to look 
at three-fold coincidences through consecutive scintillator detectors. Since energy 
loss for muons is much smaller compared to positrons in the calorimeters, muons 
can travel between consecutive detectors with little energy loss. The dashed line 
in Figure 9 is the muon losses determined from these three fold coincidences. The 
agreement between the two methods is good. The muon loss time spectrum is 
empirically added to the fitting function with a scale factor, which is a fit param¬ 
eter. 
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Fig. 8 : Average energy vs time with 
(lower) and without (upper) the pileup 
subtraction. 



Fig. 9 : Muon losses. 


It has been previously mentioned that the n was 0.137. This brought some difficul¬ 
ties to the analysis since uj c b 0 — uj a was very close to u a . The fit had a difficult time 
to separate these two frequencies. Therefore, we had a considerable systematic 
effect on a; a . However, the CBO phase changes from 0 to 2tt around the storage 
ring. For that reason, when the data from individual detectors are added, this 
effect becomes almost four times smaller. On the other hand when one looks at 
the precession frequency determined from the individual detectors, the effect is 
visible. Figure 10 shows the g-2 precession frequency obtained from the individual 
detectors. The fitting function used here was the ideal five parameter function 
including modulation of the number of detected events caused by CBO (total 8 
parameter). A fit to the data in Fig 10 gives a y 2 /DOF for a fit to a constant, 
which is unacceptable. When these data are fit to a sine wave (the CBO phase 
changes 0 to 2tt around the ring), the y 2 /DOF becomes acceptable (Fig. 10). 
The central values on u a for both fits are very close. 
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Fig. 10 : Precession frequency vs detectors when the acceptance change due to 
CBO is included in the fits. 


The next step was to make the fitting function more precise by adding the known 
effects of energy modulation. One of these effects is the modulation of the g-2 
asymmetry by CBO. The amplitude of the sine wave is reduced dramatically since 
most of the effect was removed (Fig. 11). Figure 12 represents the result when 
both asymmetry and phase modulations are included into the fit. 
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Fig. 11 : Precession frequency vs detec- Fig. 12 : Precession frequency vs de- 
tors when asymmetry modulation due to tectors when both asymmetry and phase 


CBO is included in the fit function. 


modulations are included in the fits. 















































































This study showed that the effect can be removed completely by adding all known 
effects to the fitting function. However, the center frequency value is not very sen¬ 
sitive to the type of functional form used. 


Another method was pursued in the analysis. That was to sample the time spec¬ 
trum with the CBO period, 11 so any CBO related effects can be removed from 
the data and it can be fit to a five parameter ideal function. The result of this 
method was consistent with the result of the method described in detail above. 

The described analysis determined the precession frequency with 0.7 ppm statis¬ 
tical error. The most significant contributions to the systematic error were CBO 
(0.21 ppm), pilcup (0.13 ppm), gain changes (0.13 ppm), lost muons (0.10 ppm) 
and fitting procedure (0.06 ppm). 


4 uj p Analysis 

The magnetic field B is obtained from NMR measurements of the proton res¬ 
onance frequency in water, which can be related to the free proton resonance 
frequency u p . 12 The field is continuously measured using about 150 fixed NMR 
probes distributed around the ring, in the top and bottom walls of the vacuum 
chamber. In addition to this measurement, the field inside the storage ring where 
the muons are, is mapped with a trolley device. This device carries 17 NMR 
probes on it and the measurements were repeated periodically 2-3 times a week. 

The trolley moves inside the storage ring and measures the field around the ring 
without breaking the vacuum. Figure 13 shows a two-dimensional multipole ex¬ 
pansion of the field averaged over azimuth from a trolley measurement. The total 
systematic uncertainty from the field measurements is 0.24 ppm. Figure 14 shows 
the field map in the storage ring. 
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Fig. 13 : Contour plot of multipole expansion. 



Fig. 14 : Magnetic field map. 


5 Result 

The anomaly can be obtained from the result of the independently analyzed 
frequencies a; a and u) p , and is determined to be 

„ _ UJa _ ^a/^p 

e 70 / / 

Viyl < B > [Ip/— LJa/LUp 

Figure 15 shows the comparison of the last three BNL results with the SM eval¬ 
uation. 






















































Fig. 15 : Comparison of theory and recent BNL results. 

The theoretical value of a M in the SM is determined from a M (SM) = a M (QED) + 
a^had) + a M (weak). The QED and weak contributions are given by 13 a M (QED) = 
11 658 470.57(0.29) xlO - 10 (0.25 ppm) and a /t (weak) = 15.1(04) xlO " 10 (0.03 ppm). 
The leading-order contribution from hadronic vacuum polarization contributes the 
largest uncertainty to a M (SM). Until recently, a^had, 1) = 692(6) x 10” 10 (0.6 
ppm) was the most reliable value, 14,15 where data from both hadronic r-decay and 
e + e~ annihilation were used to obtain a single value for a /t (had, 1). Recently, two 
new evaluations 16,17 using the new e + e _ results from Novosibirsk 18 have become 
available, and Ref. 16 also employs data from hadronic r-decay. While the two 
new analyses of e + e _ data agree quite well, the value of a /t (had, 1) obtained from 
r-decay does not agree with the value obtained from e + e _ data. 16 

The higher-order hadronic contributions include 19 a M (had, 2) = —10.0(0.6) x 10 _1 ° 
and the contribution from hadronic light-by-light scattering is 20 a M (had, lbl) = 
+8.6(3.2) x 10~ 10 . Using the published value of a tl (had, 1) from Ref., 14 the stan¬ 
dard model value is a M (SM) = 11 659 177(7) x 10~ 10 (0.6 ppm). 

From the most recent measurements at BNL, the muon anomalous magnetic mo¬ 
ment is determined as a fl (exp) = 11 659 204(7)(5) x 10~ 10 (0.7 ppm) 21 and the 
difference between a^(exp) and a M (SM) above is about 2.6 times the combined 













statistical and theoretical uncertainty. If the new e + e _ evaluations are used 16,1 ' 
the discrepancy is about 3 standard deviations, and using the r-analysis alone 
gives a 1.6 standard deviation discrepancy. 

In the 2001 run, the muon g-2 experiment at BNL took data with negative muons 
with a similar statistical power to the 2000 data. This measurement will provide 
a test of CPT violation and also an improved value of a /t . 
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